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Abstract We show that a VMD based theoretical input allows for a significantly improved accuracy for the 
hadronic vacuum polarization of the photon which contributes to the theoretical estimate of the muon g — 2. 
We also show that the only experimental piece of information in the r decay which cannot be accounted for 
is the accepted value for Br(r —> mrf-r), while the spectum lineshape is in agreement with expectations from 
e'*'e~ annihilations. 
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1 Introduction 

The Hidden Local Symmetry (HLS) Model[l|, Q 
implements the Vector Meson Dominance asssump- 
tion within the framework of Effective Lagrangians. 
The non-anomalous sector of this model covers anni- 
hilation channels like e+e~ tt+tt" or e+e~ — > KK 
and some important decay channels like t — + TTTrVr- 
The non-anomalous sector can be supplemented with 
an anomalous sector [il, i, i, allowin£| for 77P, 
jPV, PVV, jPPP and VPPP couphngs. There- 
fore, annihilation processes like e+e~ {n" /ri)j, or 
e+e" — > tt'^tt+tt" can enter the HLS framework as well 
as all radiative decay processes of the form V Pj 
or P^77 or also processees like 7^/77' ^ 7r+7r"7. 

Therefore, the HLS model provides a unified 
framework valid in the low energy regime up to the 
(f) mass region. It encompasses most annihilation and 
decay processes. 

However, in order to be confronted with exper- 
imental data, the HLS model should be equiped 
with symmetry breaking mechanisms. Implementing 
SU(3) breaking is done using a variant^] of the BKY 
mechanism[3| in the non-anomalous sector. Break- 
ing of the (nonet) U(3) symmetry for pseudoscalar 
mesons is also an important issue ; it is generated 
0] by determinant term Lagrangian pieces fo^. The 
SU(3) breaking of the anomalous La gra ngian is done 
following the scheme proposed by 10|, supple- 



mented with a vector field renormalization recently 
justified [li]. This full SU(3)/U(3) breaking of the 
HLS model, recalled in , has allowed a successful! 
description of all light meson radiative decays l^, M \ ■ 
A consistent treatment of the e+e~ — s- tt+tt" anni- 
hilation and the r — + tttti't decay requires an appropri- 
ate mechanism for Isospin Symmetry breaking (ISB). 
This has been defined in [12| and has improved the 
description of all processes listed above (annihilation 
and decay processes) as shown in 



2 How can VMD improve estimates of 

9-2 ? 

Therefore, the HLS model provides a framework 
able to describe in a unified way an important number 
of cross section^ with an additional set of radiative 
decay modes. These play the major role of constraints 
in order to determine numerically the parameters of 
the SU(3)/U(3)/SU(2) breaking scheme. The decay 
T — > TTTTP-r is nothing but an additional constraint, also 
subject to ISB effects usually split up into short range 
[lit ^-nd long range 19, 2C| (resp. Sew and Gbm(s)) 
corrections. These are only overall rescaling factors. 

Within this unified model all relevant data 
(already listed) depend on a very few basic parame- 
ters, namely the CKM matrix element V^d, the elec- 
tric charge e, the pion decay constant the univer- 
sal vector coupling g, the weak interaction coupling (72 
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left aside because of a still misunderstood problem concerning the ratio of two kaon decay modes of the <j) meson [l7| . 
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and a parameter named a, specific of tlie HLS modei 
0,0], expected close to 2. V^d-, fn, and 32 (related to 
the Fermi constant 32 = 2mvK Vv2Gf) are accurately 
known. Therefore, the only parameters to be fitted 
from data are a and g. The anomalous sector intro- 
duces 4 more parameters (named Ci in |2]) in such a 
way that only two parameters should be determined 
by fit [li] : the combination Ci — C2 and C3. 

The U(3)/SU(3) breaking procedure introduces 4 
breaking parameters determined by only the radia- 
tive decays 12, 1^ : Za, zy, zt and x. Some of these 
have a clear physical meaning. Indeed, za = [/k/ fn]^ 
is the squared ratio of the kaon and pion decay con- 
stants. X is the nonet symmetry breaking parameter, 
tightly related with the pseudoscalar mixing angle in 



the octet-singlet basis [l2|, 9ps — —10°. More im 



portant for the present purpose is the ISB breaking 
scheme which introduces more parameters 12, 1^ to 
be fitted and is sketched below. 

Our extended model [3| can provide a global fit 
to the whole set of data listed above. Stated other- 
wise, the parameters given above underly a physics 
content common to a very large number of annihi- 
lation and decay channels. Therefore, our overcon- 
strained paramctrization of the VMD physics allows 
for a global overconstrained fit. Then, if these con- 
straints are well accepted by the data, parameter val- 
ues and the parameter error covariance matrix will 
be defined with high accuracy. This should reflect in 
better estimates of the various contributions of the 
photon hadronic vacuum polarization (HVP) to a^. 
For this purpose, one only relies on the description 
quality of the annihilation cross sections and on the 
consistency of the various data sets with each other. 

The quality of the description of the various cross 
sections gives also a hint on the quality of the esti- 
mates these allow for a^. As stated above, the limit 
of validity of the HLS model extends to slightly above 
the 4> mass. However, this s region contributes more 
than 80% to the numerical value for and the cor- 
responding uncertainty is as large as ~ 35% of the 
total uncertainty. Therefore, even if limited, the 
expected improvements may have important conse- 
quences concerning the physics of g — 2. 



3 Breaking of the isospin symmetry: 
Vector field mixing 

Concerning the sector of neutral vector mesons, 
at leading (tree) order the ideal fields p'j, uji and 
which enter the HLS Lagrangian ~ as any VMD La- 
grangian - are mass eigenstates with resp. masses 



— m^, — and m?, — 



Zym? {rn? = ag'^ fl). How- 



-p ■■-cj ■■" ■■"4, 

ever, at one loop order, the Lagrangian piece 



iag_ r 
Aza I 



Zv<Pl 



K- dK^ 



p'j—LUi + V2Zv(t>I 



K°dK°} (1) 



induces transitions among the ideal vector meson 
fields p5, lui and through kaon loop^. Therefore, 
at one loop order, the ideal fields are no longer mass 
eigenstates and thus do not coincide any longer with 
the physical p", lu and (f> fields which, instead, must 
be mass eigenstates. At one loop order, the squared 
mass matrix for the field triplet (pj, uij, 0j) is 
given by Eq. (12) in [l3| and its eigensystem can be 
constructed perturbatively. One can define 3 mixing 
functions [l5| : a{s), (3{s), 7(5) which can be consid- 
ered as complex "angles" and are function of s, the 
squared momentum fiowing through the vector meson 
line. a(s), /3(s) and 7(5) describe resp. the p^—uj, 



p —4> and u! — (p mixings. These angles [15|, functions 
of the kaon loops and of the pion loop, contain sub- 
traction polynomials to be fitted using experimental 
data. The relationship between ideal and physical 
fields can be written in terms of these angles : 

a;-|-Q;(s) ^" + 7(3) cf) 
^ 0-/3(s) p'-i{s) oj ) 



V J 



(2) 



Therefore, the vector meson mixing induced by loop 
corrections being s-dependent, is a quite important 
feature. This transformation propagates to the in- 
teraction terms. For instance, the term describing 
the interaction of a pion pair with vector mesons be- 
comes : 



9vr+^:^[p«-a(.) c + /?(.) 4] (3) 

clearly exhibiting the origin of the isospin 1 part of 
the physically observed uj and (j) fields. 

The ^ — V transition term is also interesting. It 



lag . 



1 Actually, in the more complete Lagrangian, K*K* and KK* loops come in complementing the kaon loops along the same 
lines El. 
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becomes 



V2. 



zv 



■A- 



V2, 



Zv 



PM P° + fZ{s)oo-U{.s)^ 4> 



■A 



(4) 



where fZ{s) = ag fl[l + Ola{s} 
defined correction terms 



?(s),7(s))] has well 
The electromag- 



netic field is denoted by A. 

The most interesting feature here concerns the p 
meson which then gets different transition amplitudes 
to the 7 and W fields, One can, indeed, show that the 
amplitude ratio is : 



J p 

fW 

J p 



a{s) \/2zv 



3 



, {f:r=ag.n) (5) 



where the s-dependent terms represent the isospin 
part of the p° meson inherited from its loj and 0/ 
components. This makes different the interaction of 
the p° and fields with resp. the 7 and W gauge 
fields. 

Therefore, our isospin breaking scheme results 
in physical vector fields which are mixtures of def- 
inite isospin components and their exact content is 
s-dependent. 

4 Sketching the global fit to data 



.7r°7, 



The cross sections for e+e~ tt+tt", e+e" 
e'*'e" —^ 777 and e'^e^ — + tt'V+tt^ have been worked 
out in together with the expressions for the rel- 
evant set of decay partial widths. The expression for 
the T TiTiVr spectrum has been computed in 
and can also be found in 16|. The corresponding 



formulae have been implemented within a computer 
code aiming at performing a (simultaneous) global fit 
to all existing relevant data. 

All existing e+e^ annihilation data samples have 
been considered in the context of our global fit 
method. For the tt+tt" final state, this covers 
the former data sets collected in 21 1 and in 
and the more recent ones collected at Novosibirsk 
All existing data sets with 



23|, l2J, l25|, l26|, l27|. 



the (7r°/?7]7 final states have also been considered 



28, 29 



30, 31, 3 



For the tt^tt+tt annihilation channel, the main 
available data sets have been provided by CMD-2 
2i,[3i,[3l[3i,Ei and SND [Hi, [HI. These have been 
considered along with older data sets 38, i^l ; only 
the very old data set from [io) has been eliminated 



because it was not clear how to account precisely for 
its systematics. 

Actually, after analyzing the scale uncertainties 
claimed for the CMD-2 and SND three pion data sets, 
we were led to leave aside [15[ also the SND data sets 

3i,E3. 

Finally, the tt+tt" KLOE data set, collected at 
Frascati using the ISR method and reanalyzed re- 
cently [4l|, has been included in the data sets con- 
sidered. 

Concerning the t -kttv^ spectra, we considered 
those from CLEO 0, ALEPH and BELLE [44|. 
These data sets will be commented with some details 
in Section 6. 

Full information about the fit properties and qual- 
ities can be found in [lil . [l6l | and are not presented 
here because of lack of place. Let us only mention 
that they are always affected by very good probabil- 
ities. 



5 Improved estimate of the photon 
HVP 

In order to estimate the various contributions of 
the photon HVP to for s < 1 GeV, we followed a 
specific procedure : 

• Use always all the e+e~ annihilation data sam- 
ples essentially collected at Novosibirsk |2ll 



23 



24, 


^ 


26^27L 28 


. I29I 


30, 


31 


33, 


34, 


35, MM 









Use always the various partial widths of types 
VP^ and P — !■ 77 as reported in the Review 
of Particle Properties [lil . These play a crucial 
role in order to overconstrain our model param- 
eter values. 

Examine the effect of the tt+tt" KLOE data 



41| separately, because the fit properties of this 



sample are not fully satisfactory. 

• Add as further constraints, separately or al- 
together, the T data from C (CLEO B 
(BELLE 44]) and/or A (ALEPH 0), in or- 
der to exhibit the specific infiuence of each. 

In the comparison with experimental data, we fo- 
cus in the following on the contribution of the pion 
loop only {i.e. a^(7r7r)), integrated between ^/s = 
0.630 GeV and = 0.958 GeV. Indeed, most ex- 
perimental groups have published their estimates for 
a^(7r7r) in this reference energy range. As these exper- 
imental results are corrected for final state radiation 
(FSR) effects, we do alike. 



No. X 



M. Benayoun et al: Can VMD improve the estimate of the muon g — 2? 



4 



In order to check our method and iUustrate its ef- 
fect, we have first run our code using each of the data 
sets from 2^, 25| and 27| in isolation, together with 



our full set of radiative decay information (17 pieces). 
The results derived from the fitted pion form factor 
are reported in the first 3 lines of Table [T] and the 
errors shown are the total uncertainties. Indeed, the 
fit is done with a procedure combining appropriately 
ISl . Il6j statistical and systematic errors. 



Table 1. Our estimates for 10^" a,i(7r7r) and the 
corresponding experimental values from resp. 
[i^] . [iit, [irt - The last two lines are averages 
proposed by [iit . 



Data Set 


Exp. 


Value 


Reco. Value 


Prob. 


CMD2 (1995) 


362. 1± 


2.4±2.2 


362.9 til 


51% 


CMD2 (1998) 


361. 5± 


1.7±2.9 


362.2±2.1 


49% 


SND (1998) 


361.0± 


1.2±4.7 


361.0±2.1 


99% 


NSK (all) 


360.2± 


S.Otot 


361.7±1.3 


48% 


NSK +KLOE 


358. 5± 


2Atot 


362.1±1.1 





One clearly observes an important improvement 
of the accuracy following from having built, for the 
first time, a working model which simultaneously fits 
the radiative decays and the annihilation data. Com- 
paring the results obtained using each of the CMD-2 
and SND data sets in isolation and altogether, the ef- 
fect expected from an increased statistics is observed 
with its expected magnitude. The net effect is a fac- 
tor of ~ 2 improvement of the uncertainty. As will 
be seen shortly, this is also due to the fact that the 
uncertainties (and biases) within the data sets just 
quoted are well under control. The last data column 
gives the probability of the underlying fit to the pion 
form factor and the decay data. The fit probability 
of the SND data clearly refiects a too conservative 
estimate of their systematics. 

KLOE data [41:] help in slightly improving esti- 
mates at the expense, however, of a poor fit proba- 
bility, essentially due to a (still) poor control of the 
systematics within this data set. 

Our favorite estimate of (tttt) (fourth line in Ta- 
ble [T]) compares favorably to the newly issued experi- 
mental results produced from recent pion form factor 
data collected using the ISR method by the KLOE 
Collaboration (lO^o a^(7r7r) = 356.7± 0.4± 3.1) 
and by the BaBar Collaboration S| {W° a^^Tn:) = 



365. 2± 2.7). These two new measurements illustrate 
that one needs motivated theoretical input in order 
to take a full profit of the new high statistics data 



sets. Indeed, [48[ proposes an average of the four ex- 
perimental values given in Table [Tl fjii. [jj. [27ll41|) 
and of the BaBar estimate over the same energy 
range ; using a sophisticated statistical method, they 
get 10^° a^(7r7r) = 360.8 ±2.0tot. 

Comparing this average with our fit value (last 
line in Table[T]) - which does not use the (not yet pub- 
lic) BaBar data - is interesting. Indeed, it shows that 
the increased statistics provided by the ISR method 
at DAPHNE and BaBar has not allowed a real break- 
through in the accuracy of (tttt) , because of the sys- 
tematics specific to each experiment and of the diffi- 
culties encountered while merging the different data 
samples. 

Instead, what is illustrated by Table [T] is that an 
adequate theoretical input - like VMD - may allow 
sizable improvements. Of course, the relevance of 
this input should be (and actually is) reflected by the 
global fit qualities 



15 



We do not discuss here the effects of introducing 
the {tt° /r])'y and Tr+TT"?!-" data ; this has been analyzed 
in full details in 



16| 



Let us, nevertheless, mention 
that these data sets, with poorly known sytematics, 
allow to confirm the central values for (tttt) without 
a visible improvement of its uncertainty. 

6 Adding the r spectra to the fitted 
data samples 

As mentioned in Section 4, we only deal with the 
CLEO (C) 0, ALEPH (A) and BELLE (B) 
4J] data sets. The (C) data set provided by CLEO 
is actually the normalized spectrum 1/NdN/ds. The 
absolute normalization for dT{T ■n'Rv^) / ds is deter- 
mined by a multiplicative factoiH, the branching ratio 
Br(r rnrVr)- Therefore, the CLEO spectrum we use 
is not sensitive to this branching ratio. As, following 
the BELLE Collaboration 4j], we allow for a rescal- 



ing of the B data set, we are only marginally sensitive 
to Br(T — + tttti't). Instead, as there is no reported 
uncertainty on the normalization of the ALEPH (A) 
|i^.^(s)P spectrum, we have not allowed any rescaling 
for the A data set. 

This way to proceed with B and C is not the usual 
one. Indeed, usually, the B and C |i^^(s)P spectra are 
constructed as their reported normalized spectrum 
1 /NdN/ ds multiplied by the world average valu«0 for 
Br(r tttti/t-) 



44 



In the (global) HLS model, the r spectrum is de- 



§See, for instance, Eq. (7) in |4J 
IT As can be concluded from Figure 6 in [46 
corresponding ALEPH |43| measurement. 



the world average value for Br(T — > ixTrVr) differs only marginally from the 
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termined essentially by the Higgs-Kibble mass 
(occuring in the Lagrangian) and by the cou- 
pling to a pion pair. Naming the p° mass squared 
m^{— ag^ fl) and g its coupling constant to a pion 
pair, we have defined the corresponding quantities for 
the meson by m'^ + Sm^ and g + Sg. Interestingly, 
the absolute magnitude of the r spectrum and the 
yO* width are both determined by the p'^^tttt coupling 
constant and then by g+Sg. Isospin symmetry break- 
ing effects specific of the r decay modify this picture 
by introducing short range [l^ (Sew) a-nd long range 
UM (G'bm(s)) corrections which both factor out and, 
therefore, contribute to the absolute magnitude of the 
T spectrum. 



CMD2 (1995) 
CMD2 (1998) 
SND (1998) 



NSK new timelike 



NSK all timelike 
+ n° /ri 7 

+ +7t* 7T" tt" 

NSK+B+C 

NSK + B'+C 

NSK+KLOE 

NSK+KLOE+B+C 

NSK+A 

NSK+B+C+A 

NSK+KLOE+A 



345 



350 



a, (n*,n")*10'° 



Fig. 1. Figure 1. The tttt loop contribution to 
g — 2 integrated between 0.630 and 0.958 GeV 
using our global fit running with various com- 
binations of data sets. 

We have first performed fits with e~^e~ and r data 
in order to determine 6m^ and Sg. It happens 1^ 16 1 
that the fits return Sm^ and dg consistent with zero 
at a ~ 1 (T level. Therefore, we do not find signifi- 
cant differences between the p'' and (Lagrangian) 
masses and couplings. It thus follows that the differ- 
ence between the pion form factor in r decays {F^ (s)) 
and the 1=1 part of the pion form factor in e'^e~ 
annihilations {F^{s)) is fully carried by the factor 
a/ SewGem{s), which affects the r dipion spectrum. 

Then, fixing 5w? = Sg = 0, we have redone our fi- 
nal fits allowing for a rescaling of the B data sample. 



by varying the set of data sets (listed in Section 4) 
submitted to the global fit. 

7 A localized failure of CVC ? 



Our global fits are always fairly good 1^, 3| and 
result in an overall rescaling factor for the B data 
sample 1 + A with A= (-4.84to 92)%, in good corre- 
spondence with the BELLE fit result 4j| which can 
be written Xbelle = — (2±1±4)%. In this approach, 
the C and B data samples are always well described ; 
the ALEPH spectrum is reasonably well described be- 
low 1 GeV, however more poorly than the C and B 
data samples [l^. This is partly du^ to the fixed 
absolute normalization of ALEPH data, i.e. to the 
accepted value for Br(T ^ tttti/^). 

At this step, one should note that the HLS model 
we use, equiped with symmetry breaking schemes ac- 
counts fairly well for : 

• all e'^e~ annihilation cross sections listed in Sec- 
tion 4, 

• all partial width decays of the form P77, VPj 
and 77/77' —> 7r+7r~7, 

• the lineshape of the dipion spectrum in t decay, 
especially those provided by CLEO and BELLE 
which are quite similar. 

Stated otherwise, the single piece of information 
which does not fit within this overall picture is the 
accepted absolute normalization of the t dipion spec- 
trum, i.e. Br(r — > 7r7rz^.r)- 

If one excludes an experimental bias, one thus 
needs a specific additional breaking effect affecting 
solely the t decay. However, as this missing piece re- 
sembles a global rescaling of the r spectrum, a possi- 
ble candidate could be a revised Sew factor numerical 
value, if relevanlF*!. 

8 Influence of the r spectrum 

Figure [T] displays the value for a^(7r7r) integrated 
along the canonical interval around the p peak, as 
coming from our (global) fits. The 4 upmost data 
points are the values shown in Table [TJ Thus, the 
fourth line gives the result derived from a combined 
fit to the data given in The fifth line 

displays the result coming from the combined fit to 



" Indeed, if the CLEO and BELLE pion form factors are in fairly good agreement with each other, they sensitively differ from 
the ALEPH form factor in the very low and in the high energy regions, as can be seen from Figure 12 in 14411 . 

**It is generally assumed that the value found for t^ttu 18] coincides with the corresponding factor for t^ttttu. This crucial 
assumption does not seem to have been proved. 
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the 7r+7r~ data sets just quoted and to the older tt+tt" 
data sets given in [2l|, [22] ■ For the hne indicated by 
+7r''/777 , we have added the corresponding data sets 
to all 7r+7r~ data. In order to get the result indicated 
at the line flagged by ++7r+7r~7r°, the corresponding 
data samples have been considered together with all 
the previous ones. Concerning the rest of Figure [TJ 
NSK denotes all e+e~ annilhilation data combined 
with KLOE, ALEPH, BELLE, CLEO in the way in- 
dicated at the corresponding line. 

One can conclude from Figure [1] that all data set 
combinations submitted to fit and built up from all 
e+e~ annihilation samples and from the B and C sets 
provide quite consistent results. Instead, as shown by 
the 3 downmost a^(7r7r) values, including the ALEPH 
data set always provides a shift upwards by ~ 5 10~^". 
This is almost certainly related with the branching 
ratio issue discussed in Section 7. 



9 Conclusions 

We have proved that a theoretical VMD input 
permits to significantly improve the accuracy of pre- 
dicted value for the muon g — 2 value, as clear from 
Table [TJ Some further improvement is reached by 
adding the t spectra, however marginal. Our VMD 
input certainly increases the disagreement between 
the expected value for the muon g — 2 and its direct 
BNL measurement. 

Another important remark is that the p meson 
lineshape observed in the t dipion spectra in perfect 
agreement with expectations from VMD. The single 
surviving issue in our data set, the largest one ever 
analyzed within a single model, is solely the value 
for Br(r mri/r), expected slighly smaller than its 
presently accepted value. If not an experimental bias, 
this may indicate that symmetry breaking effects in 
T decays are still to be revisited. Until this issue is 
clarified, one should consider cautiously the predic- 
tions for the muon g — 2 provided by the r dipion 
spectrum, especially those depending on the absolute 
scale of this spectrum. 
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